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SUMMARY 


This  report  describes  a  computer  code  which  calculates  a  geometry 
parameter  from  inputs  of  an  air  vehicle  aft  end  area  distribution.  This 
parameter,  called  the  "Integral-Mean-Slope"  or  IMS,  was  found  to  correlate 
well  with  high  subsonic  and  transonic  aft  end  drag  measurements  and  is 
used  a_  a  predictor  of  pressure  drag. 
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INTRODUCTION  AND  BACKGROUND 


The  NAVAIRDEVCEN  was  tasked  with  an  assignment  from  the  Naval  Air 
Systems  Command  to  develop  an  accurate  method  of  predicting  aft-end 
fuselage  drag  of  twin  nozzle  "fighter"  type  aircraft.  A  computer 
program  (reference  (a))  and  several  supporting  documents  (references  (b) , 
(c)  and  (d))  were  developed,  in  addition  to  this  report,  as  a  result  of 
the  aft-end  drag  analysis  methods.  During  the  methodology  development, 
it  was  clear  that  one  of  the  major  parameters  which  correlates  well  with 
the  aft-end  drag  was  the  Integral  Mean  Slope  (IMS) . 

This  report  presents  a  computer  code  which  allows  the  user  to 
compute  the  IMS  from  inputs  of  the  vehicle  geometry.  The  input  format 
is  compatible  with  other  in-house  graphics  utility  codes  such  as  refer¬ 
ence  (g) ,  so  the  user  may  preview  the  inputs  before  execution. 

The  IMS  parameter  is  calculated  by  integration  of  the  normalized 
rate  of  change  of  fuselage  area  as  defined  by  the  following  equation. 


1.0  r 

A  /A  J  . 
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where. 


Ag  represents  the  aft  area  integration  limit 
A^^  represents  the  maximum  fuselage  area 
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D„„  represents  the  maximum  equivalent  fuselage  diameter. 
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The  IMS  parameter  was  originally  introduced  by  Pratt  &  Whitney  Aircraft 
in  reference  (e)  as  a  correlation  parameter  relating  fuselage  area 
distribution  with  the  variation  of  aftbody  drag  coefficient.  Another 
variant  of  the  IMS  is  the  truncated  IMS  or  IMST  which  was  formulated  by 
Boeing  in  reference  (f).  The  truncated  IMS  parameter  is  a  recent  varia¬ 
tion  of  the  basic  IMS  calculation  which  utilizes  empirical  data  associated 
with  probable  streamline  separation  '■s  a  function  of  aftbody  curvature 
and  Mach  number.  The  calculation  procedure  involves  testing  the  local 
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aftbody  slope  at  a  given  free  stream  Mach  condition  for  the  possibility 
of  separation  and  the  modification  of  the  aftbody  slope  to  conform  to 
the  separation  streamline  when  appropriate.  The  empirical  data  asso¬ 
ciated  with  the  truncated  IMS  parameter  is  presented  in  Figure  1  and  was 
extracted  from  reference  (f). 

CODE  OPERATION 

The  code  consists  of  a  main  program,  IMS,  and  two  subroutines, 

SPLNQl  and  SL0PE1.  IMS  processes  the  input  data  to  create  an  array  of 
the  normalized  area  distribution  versus  normalized  aftbody  length  (see 
equation  (1)).  An  integration  is  performed  in  accordance  with  equation  (1) 
above  in  steps  of  .01  of  the  normalized  aftbody  length  to  arrive  at  the 
IMS  and  IMST. 

The  subroutines  SPLNQl  and  SL0PE1  associated  with  this  computer 
code  utilize  a  spline  fit  curve  through  the  input  data  to  produce  a 
piecewise  cubic  with  continuous  first  and  second  derivatives.  The 
aftbody  area  distribution  is  curve  fitted  using  the  SPLNQl  subroutine 
for  interpolation  of  data  values  and  the  calculated  aftbody  slope  is 
obtained  in  a  similar  manner  using  the  SL0PE1  subroutine. 

OPERATING  INSTRUCTIONS 


The  IMS  parameter  correlates  well  with  high  subsonic  and  transonic 
aftend  drag  measurements  (see  references  (a)  and  (b))  if  sufficient  care 
is  exercised  in  the  precision  of  input  data.  The  numerical  calculation 
procedure  necessitates  a  smooth  input  data  curve  of  aftend  area  dis¬ 
tribution  for  maximum  precision  of  calculated  IMS  parameter.  Therefore, 
it  is  advisable  to  plot  the  area  distribution  prior  to  use  of  this 
program  to  insure  a  smooth  input  curve.  The  data  input  format  is  com¬ 
patible  with  the  NAVA1EDEVCEN  plotting  computer  codes  (reference  (g))  to 
facilitate  the  input  of  analytically  smooth  data  curves. 

The  program  input  data  requirements  are  specified  in  Table  I,  where 
a  minimum  of  nine  data  cards  are  required  for  execution.  The  maximum 
integration  limit  (XMAX)  is  input  in  the  first  ten  columns  of  card  1  and 
can  represent  either  fuselage  station  directly  or  can  be  normalized  with 
the  maximum  equivalent  fuselage  diameter.  The  number  of  Mach  numbers 
(XNOMACH)  is  input  in  the  next  ten  columns  and  the  maximum  fuselage  area 
(AMAX)  is  entered  in  columns  21  thru  30. 

The  fuselage  area  distribution  (input  on  card  7)  can  be  specified 
in  dimensional  units  of  either  square  feet  or  square  inches  by  setting 
the  parameter  (ADIMEN)  in  columns  31  thru  40  to  either  0.  or  1.,  respec¬ 
tively.  The  user  may  also  input  the  normalized  area  distribution  as  a 
function  of  normalized  aftbody  length  X/DEC,  provided  Amax  is  input  as  0. 
Detailed  integration  steps  during  the  process  of  calculation  of  the 
overall  IMS  parameter  can  be  printed  by  insertion  of  a  value  of  1.  in 
columns  41  thru  50. 


--r.-.as.. 
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FIGURE  1  AFTBOOY  SEPARATION  CRITERIA 
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INPUT  DATA  FORMAT  FOR  IMS  PROGRAM 


OE-IZ 

samvA  viva  xvwa/x  oz-n 


TABLE  I  (CONT*D) 

INPUT  DATA  FORMAT  FOR  IMS  PROGRAM 


BLANK  CARD 


NADC-76 386-30 


The  second  data  card  is  reserved  for  Input  of  the  Mach  number 
values  (XM)  as  real  numbers  in  fields  of  10  columns  starting  in  column  1. 
The  third  data  card  allows  the  input  of  a  case  identification  or  title 
(IDENT)  with  cards  four  and  five  input  with  blank  fields.  The  number  of 
fuselage  data  stations  (NFP)  is  input  on  card  6  as  an  integer  value 
right-adjusted  in  columns  5  thru  7  with  the  fuselage  data  values  (FX) 
input  as  monatomically  increasing  real  numbers  in  fields  of  ten,  starting 
in  column  11.  If  necessary,  additional  fuselage  data  values  can  be 
continued  on  subsequent  cards,  starting  each  in  column  11.  The  number 
of  fuselage  area  data  values  (NPP)  is  input  on  card  seven  as  a  right- 
adjusted  integer  in  columns  5  thru  7,  and  must  correspond  with  the  value 
input  on  card  6.  The  fuselage  area  data  values  (FX)  are  input  in 
correspondence  with  the  previously  input  fuselage  values  as  real  numbers 
in  fields  of  ten  starting  in  column  11.  Additional  fuselage  data  values 
can  be  entered  on  subsequent  cards,  starting  each  in  column  11.  The 
final  card  for  each  case  is  input  with  blank  fields  and  additional  cases 
can  be  Included  by  repeating  the  input  format  of  cards  1  through  8 
followed  by  a  blank  stop  card. 

The  output  provided  by  the  program  consists  of  the  standard  IMS 
parameter  calculation  and  the  truncated  IMS  parameter  for  each  Mach 
number  specified  by  the  user. 

A  sample  case  is  included  in  Appendix  A,  together  with  the  re¬ 
sulting  sample  output.  A  listing  of  the  IMS  source  program  and  required 
subroutines  are  included  in  Appendix  B. 
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In  this  sample  problem,  the  Inputs  on  page  A-5  show  that  the  aft 
body  area  distribution  was  input  in  normalized  form  using  17  data  sets 
terminating  at  a  normalized  aft  body  length  of  2.611  equivalent  diameters. 
From  the  inputs,  we  can  see  that  the  IMS  was  to  be  determined  at  5  Mach 
numbers  (.5,  .6,  .7,  .8,  and  .9)  from  the  maximum  area  station  +.  a 
length  of  2.611  diameters.  The  IMS  computation  can  be  terminated  at  a 
more  forward  body  station  by  specifying  the  parameter  XMAX  to  be  less 
than  2.611  in  this  case. 

The  output,  on  page  A-6  shows  the  value  of  IMS  and  IMST  to  be 
.63665  over  the  Mach  number  range  indicating  no  probable  separation  by 
the  criterion  discussed  in  the  introduction  section. 
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APPENDIX  B 


COMPUTER  LISTING 


The  listing  provided  herein  includes  the  IMS  program  with  subroutines 
SPLNQ1  and  SLOPE1  which  are  currently  operational  on  the  NAVAIRDEVCEN 
CDC  6600  computer  system. 
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IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  004 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 


RCAO  100. RPP.  (FIT?  J  .  I  * W  .  U  £  I 
J  •  4t  *  ' 

I  r I RPP  GT  7) RCAO  t  10  .  I  FI  (  I  i  .  I  > 

I C  A  3  10 .  IOUM 
F |  i IP*  I  ’  •  0 
F I  1 RP*  J  J  »  - 

RP  *  RPP  •  1 

pair*  iso.mii 

PRINT  ,20  ,  t9Ui.i»9.»»> 

FORMAT i  iO* . 7F 10 .4  I 

ISTRT  .  f«i2' 

IMAll  •  ; IMAX -ISTRT t»OF AC 

1 1 1  me  eg  < .  i  co  to  iao 
DC  '3D  1 <  2  . RP 

lllll  •  ( F| ( I > -ISTRT) «0F AC 


RP 


RP  • 
RP  • 


RPP 

F| | I)  . 1 *4 . RP  I 


PRINT  120 
PRINT  AO 
RRIRT  ISO 

FQRMAMTi  IPMACH  rO«  .  T224TRURCPTCD  IMSa.TAr.STO  !**•! 
IFIAFAC. 10.1.)  GO  TO  » ?0 

oo  no  in.oo 

mi)  aafac 

r  *  r . «m 


F*  1  1  I 
00  770 
RMS 
IMS 

93  200 
1*0 
»"Si 
OlSOi* 
AREA 
T  1 

<s 


*  I . C20 


SPlKOt ( l . F|. l) 

r  i 


siom  i  i  ,r«.  *.»ouT2J 

l  10  I  F  (  (  1*0  1  1ST  1M*H>  01*  |M|  |  X  •  I 
I  •  1*01 

T2  ■  SPCMOl (  I , Fl , I ) 

TPJ  •  SL0*C ifi.FI.  I . V  OU  T  2 1 

*•(*•1.1? . IMS)  VP2*1MS 

ARIA  •  ARrA*(r»i*rR7I*. *4Tf7'Pf » 

J F t IRRIRT . R( .0. |  PRINT  ’ *0 . ▼ 2 . * • 24 A t A 
1*0  lORMAtl*  AAMAI<<C<2  ».#  SkOPC>#.C<3 

'  *  r  * 


•  ARIA*#  ,  C  1  2  .  S  ) 


*P2 

;r  Bi  M  ot»>  dc  to  no 

RM*  .  ARC  A 

JOO  IMS  •  *SPl R«  t  I  •  , 2  - 1M1  I !  ' 

RRI R*  2 1 0 . IM ( I  1  . RMS. RM*. 

110  FORMA* l *1 . F i o  R.Sl.FtO  *.**.*'*  S' 

3>D  rORTJRul 
•RlRt  230 

2*0  FORMAT | 11 . v  r///| 

IF(RO»T .*«. SI  0001*0 

GO  ?0  10 

t»0 

FUNCTION  SH ■* i  (RIOC . t . IIROCOl 
C  CUll  C  Ilk  INI  »M  RIVIIIO  10/21  ”  M  CAOO* 

C  T  #•  I  S  FfRSJOR  HA*  0*0  0*T1D»  MMIRI  ASi  OF  TMC  SPktRC  COCFF|CItRTS 

C  ARC  COMPUTE*  ARC  S?0*tO  |R  TM|  ARRAY.  FOR  R  ORTA  POINTS  3#R*J 


IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  00 I 
IMS  001 
IMS  001 
IMS  00< 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  00* 
IMS  007 
IMS  007 
IMS  OOT 
l M S  007 
IMS  007 
IMS  007 
IMS  007 
IMS  OOT 
IMS  907 
IMS  OOT 
IMS  00* 
IMS  00* 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOP 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
IMS  OOS 
(M$  OOS 
IMS  OOS 
IMS  OOP 
IMS  OOS 
(MS  OiO 
IMS  OiO 
IMS  OiO 
IMS  0 10 
IMS  0 \C 
IMS  OiO 
1  **S  0-0 
l MS  0»0 
:m*  o<o 

IMS  C  D 
IM*  Oi ’ 
IMS  Oil 
IMS  on 
t ms  on 

S*  k  ROOO 
SRkROOO 
SPiROOO 
5*1 ROOO 
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STORAGE  LOCATIONS  Ai|  REQUIRED  'OR  Tm| 

ORTA  ARO  THE  COEFFICIENTS 

SNINOOO 

Mix  MATURE  IS  QUICK  lOOk-ua  fO*  LARGE 

ARRAYS 

SNL  NOOO 

dimension  cmoo i .  sot  too ; .  < <  1 1 

SNl ftOOO 

*1 R*4I aOEN 

SNl  NOOO 

"S*»lOC 

SNINOOO 

ROOT  $• I ( NS  J 

SNl NOO  i 

; a>«s*«oAT  s 

SNl ROO  ' 

mSN  t • A*« ■ 

SNl ROO  I 

NS»2*RS*2 

SNL  NOO  1 

hS2*mONT  s* 1*«SR i 

SNL  NOO  - 

1 » 1 1 »S3  t 

SNl NOO  > 

i SC • AS3* 1 

SNV  nOO  i 

1 RHODE* 1 ( L  SC  1 

SNl NOO  1 

R  ■  l 

SNl NOO  1 

•  l  •*»»' 

SNl  NOO  I 

MM*  10 

SNl NOO  7 

NT  RAO* • 1 

SNl ROO 3 

I' lN00ts.ll. n  GO  TO  130 

SNl ROO  3 

IMIIMI  10)130.  »30 . 10 

SNl NOO  3 

10 

NT  RAO*0 

SNIR002 

GO  to  <10 

SNl ROO 3 

JO 

IMUR-ll  ASN*  1  130.40. 40 

SNl NOO  2 

10 

NT  RAN* 1 

SNl ROO  2 

40 

«-*IN2 

SNl  ROO  2 

GO  TO  140 

SNiNOoa 

SO 

If  ID  no. « to.  to 

SNl  NOO 3 

40 

2  f  t  It  ••MR  I  170.140.10 

SNl  NOO  3 

70 

RM*K 

SNl N003 

K  •  K  *  1 

SNl NOO  3 

•0 

miiN-mi )  too.  no. to 

SNLNOOS 

to 

Rl  •« 

SNl R003 

GO  TO  MO 

SNIN003 

too 

NH  •  R 

SNtROOJ 

t  10 

«• tNN-Nl ) / 3*  Nl 

SNl R003 

I F 1 R • Ni ) 10 . 130.40 

SNL  N003 

1  30 

7 OUT  *  4 ( A  SN  3  ) 

SNINOOA 

CO  TO  J40 

SNl N004 

1  10 

N  *  NN 

SNl N004 

140 

NU 

SNIH004 

41  III! .» 

9NLH004 

M*X*mONT  S 

SNl R004 

If  U«  IQNODE  .  C T  .01  GO  TO  310 

SNl ROO  4 

t  SO 

I  2  ■  R I R  SR  <  ) 

SNl H004 

13*4 ( RSR31 

SNl R  0©  4 

133*43*13 

SNl  nOO • 

*3*41  10*21 

$#l NOOS 

»3243 2-  1 YJ - I ( IO*i ) 1/432 

SNl NOOS 

6111*0. 

SNLNOOS 

sim •-  s 

SNLNOOS 

Rl ■ RON  T  S  *  1 

SNlROOS 

!f(Rt  .IT  . 31  CO  TO  170 

SNiHOOS 

DO  140  1*3. N1 

SNl NOOS 

«•  *  NSN  1  •  I 

SNl NOOS 

SNLNOOS 

Il*R3 

SNlROOS 

43*13 

SNlROOS 

431*433 

SNlROOS 

43*4(41 

SNl NOOS 

133*13*12 

SNlROOS 

73*73 

SNl NOOS 

*  3  •  4  (  K  1  I 

SNl ROO  A 

*21 13 l *733132 

SNl NOO  A 

73  343 3*  173  **3) /133 

SNl N004 

x*i  . /t 2  • (  43  -  1  1 1 -431**11! * i 1  1 

SNL  R004 

Sllli * I 7  ? • w 

SNL  MOOS 

1  40 

G  !  1  )  •  W*  (  4  •1732132*72H3,I*R2I*C(1*1M 

SNl MOOT 

1  to 

CN1  *GI  R*  WT  2  •  St  •  >  1  1  > 

SNl ROOT 

If  U -CT  0  1  GO  TO  <10 

SNl ROOT 

10  1  • NO N  T  S 

SNl ROOT 

KRAS*  ROOT  S*l SC 

SNl ROOT 

1  IKRAS1 *KH1 

SNl ROOT 

GO  TO  ifO 

SNl ROOT 

1  ao 

5  0 1 ■ l 0*  2 *X 

SNi ROO  T 

no 

00  300  1*3. tot 

sniroot 

CM3*CMt 

SNl ROO  T 

INI  *G | Nil  - Sll N  t  1*1X3 

SNl 4004 

UNfiSCI  *1X1 

SNl NOOI 

300 

N  1  •  N 1  •  t 

SNl NOOI 

If U0T.01  GO  TO  330 

SNl NOO 1 

3*0 

RSX*NS3*X-NS*I 

SNINOOI 

KKi *4 l*SH- t  t 

SNl NOO • 

CXJtlIRSX) 

SNl NOOI 

330 

$•« (Ml *I(X> 1 1 

SNl NOOI 

If ( NT  RRN  IT  91  CO  TO  230 

SNl NOOI 

II«X-RTRA» 

SNl BOO* 

I 7  *  1 « ♦ rONT S 

SNl M004 

ISM!  N 

SNl MOOt 

I  I  N  ■  4  :  1  1 

SNl MOOt 

I  1  M  (Mi  •  1  {  H 

SNlROOt 

22*11  RM1X-  i  , 

5  # l MOOt 

70UT. t  f  (lX2*22*I3  fX1*I1*1 1  !•  5*  4  '  A  1  * 

4  1  A  -  *  1  </ $ 

|Nl NOOt 

1  - l€XJ-|Xt I4S/S. I#  t «$'4INt*l 1  17! 

SNlROOt 

GO  TO  340 

SNlROOt 

330 

23*I1R*41X-1 1 

SNl NOOt 

YQUt • m 1X2  *  IX  f )#22*3  *S*£X» ) • 2  3 • 1  •(> IN) *R(R*  t 1 ) -S*S* 1 IXl*|Xi* 

SNlROOt 

1  CX2I >#28/(4  *S1*1(N* M 

SNl R010 

3«0 

SNl HR i • 7  OU T 

SNl NO  T 0 

RETURN 

$Nl NO  1  0 

f  NO 

SNlROtO 

FUNCTION  Si  ONI  1  ( Rl CC . I . 1 ! N0| N , TOUT  3  1 

cuite  inure  fit  revised  <©/3m,  «  caoot 

000 

THIS  VERSION  HAS  RaO  ONTIOR  XMflt  III 

OF  THE  INURE  COEFFICIENTS 

000 

ARC  COXNuTtD  ANO  STORED  !R  TMf  ARRAY 

FOR  N  OATA  NOINTS  3»R*3 

000 

STORAGE  locations  arc  required  TOR  the 

OATA  ARC  •*£  COEFFICIENTS 

000 

REX  FEATURE  IS  RUlCR  iOOK*UR  FOR  LARGE 

ARRAYS 

000 

OIXERSIOR  Gl  >00'  S»i  <001  .11  i  1 

000 

IIR*t|ROE* 

000 

NS'RiOC 

000 

RONTS*! 1 Rl 1 

OOi 

I 0 • R  S*  NO  NTS 

00  < 

■  S#  '  *  NS*  i 

00' 

•SN2*RS*2 

OOi 

00  1 

l • R ( RS2 1 

00* 

LSC*NS3*i 

001 

I  RHODE • 1 1 1 SC  1 

00  l 

4*1 

001 

NL  *N|NI 

00  1 

NH  *10 

003 

Nt ran*  * t 

00} 

o 
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in  moms*  iMie.ni.it 
to  if  1 1 1  i  ( i d 1 1 10 . it* . ao 

20  l'lif'0 

60  10  <40 

)0  MO  .  »0 . 10 

40  irilfM 

to  «*»in 

60  10  Ito 

•  0  1* ll J 1 ?«.  I  30 . >0 

to  1 1 1 1 1«- *«* n to . ! so .  ioo 

•o  ■ 

*»*•  i 

so  inui'iiiiMte.iie.ieo 

100  ll*ft 

CO  10  12© 
f  10  »*•« 

)  26 

It t*-»l  t 90 .  >*0. *0 
ll«  r0Uf*l<»i*2) 

CO  10  1 SO 

1 40 
HO 

!(■*?»•* 

«*H«40MS 

If  U*]Q«OOC  I  HO.  1  *0. 2)0 

Ito  1 2 • 1 1 ftt* 1  I 

nmistit 

IlMl  tJ 

*i«iiio«ii 

61*1*0 

tinn-.j 

l"<OMl  i 
OO  i  TO  I'M' 

4  *  OSP 1*1 
01 •J*«0»lt 
IMII 
1M) 

III >132 
II'IUI 
m*i2>i2 

TJ.VJ 

f Ji *912 
rjfrjrt 

■  •HMD/]  -i2i«**U*  il /«• 

Sf U i ■ l 12/(*«t  I 

110  6II»»l*J2/lJ2*T2t/«2>*12«*C(I-H/t  l/K 

(«'  ■G(f  1/  12  ♦$*<■'  I  > 

iruiiie.iit.no 

190  I  0 ’  -  HO  H  T  S 

I«ts«l0tis*lic 

KRQttl  •  t  H  i 
CO  10  200 
190  l0»MO*2-« 

200  00  2 10  1*2. 10i 

CM2-CH1 

CHi »6(0l !**•<■> /  *<H2 
iiainiCMtM' 
no 


oo? 
002 
00  2 
00  2 
002 
oo? 
00  2 
00? 
00] 
00? 
003 
00] 
00] 
003 
003 
003 
003 
003 
004 
004 
004 
004 
004 
004 
004 
004 
004 
004 
OOS 
000 
OOS 
OOS 
OOS 
eos 
OOS 
OOS 
OOS 
OOS 

oot 

oot 

oot 

oot 

oot 

oot 

oot 

oot 

oot 

oot 

001 

001 

001 

001 

001 

001 

00? 

007 

oot 

oot 

oot 

oot 


If U 1220.220.230 
270  4S»*»S7*H*  **• i 

(Hi *11 MSM- l I 
(H2«IINSH) 

230  »■«!*• I « 

If  I HTttt 1 2S0. 240 . 24© 

240  ll*H-ftf*44 
I  »  ■  I  1*11007  S 

Zl'«(Kl *»fft 

» i 

VOUt • ( i ( (H2»Z2*Z2  *  fH • • 1 1 • I ! 1/  2 .  ♦  *<ft|-*C«-ll</S 
I  - l (H2 • ( H i i*S/t .  ) 

T0W?2  •  0 

60  ro  2 to 

2S0  U'llfUtDl 
2 l ‘1 (H) - 1! ft 

VQtlt  •  ( t  ( IH3*  lt4I2*IIH*l*olr  I  /  2  -  ♦  I  )/S 

1*1  f>!2  *  tH  I  I  •  S/ 1 .  > 

TOUT  2  •  I  (Nt«l3*C#M4*«  \t\ 

2 to  SiOO|i*r0UT 

miuM 

(NO 


oot 

oo* 

oot 

oot 

oot 

oot 

oot 

004 

00* 

009 

oot 

00* 

oot 

00* 

00* 

00* 

oot 


0  10 
0  10 


r-;7 

jw-iall 


fiTJC  i^-n?  not 
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